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Abstract. 3D MHD simulation of accretion onto neutron stars have shown in the last few years that 
the footprint (hotspot) of the accretion flow changes with time. Two different kinds of accretion, 
namely the funnel flow and the equatorial accretion produced by instabilities at the inner disk, 
produce different kinds of motion of the hotspot. The funnel flow produces hotspots that move 
around the magnetic pole, while instabilities produce other hotspots that appear randomly and move 
along the equator or slightly above. The angular velocities of the two hotspots are different, the 
equatorial one being higher and both close to the Keplerian velocity in the inner region. Modeling 
of the lightcurves of these hotspots with Monte Carlo simulations show that the signatures produced 
in power specra by them, if observed, are QPOs plus low frequency components. Their frequencies, 
general behavior and features describe correctly most of the properties of kHz QPOs, if we assume 
the funnel flow hotspots as the origin of the lower kHz QPO and instabilities as the origin of the 
upper kHz QPO. 
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THE KHZ QPO PHENOMENON 

The lightcurves of X-ray binaries often show a number of variability phenomena, that 
populate the power spectra of their lightcurves with different kinds of noise and oscilla- 
tory components. They are normally fitted with Lorentzian curves [1] and classified in 
terms of their characteristic frequency v and their width A v. Depending on the quality 
factor Q = v/Av, features are classified as broadband noise (Q < 2), quasi-periodic os- 
cillations (QPO, Q > 2), coherent pulsations (Av smaller than the frequency bin). Many 
weakly-magnetized accreting neutron stars show a particular kind of QPOs, whose fre- 
quencies range from 300 to 1300 Hz, called for this reason kHz QPOs [2]. Very often 
kHz QPOs appear in pairs, with the two peaks, labelled "upper" and "lower", at frequen- 
cies differing by about 300 Hz. Lower and upper peak generally differ for their proper- 
ties, the lower peak being normally more coherent and appearing in a shorter frequency 
range [3, 4]. Their frequencies are correlated with the ones of lower- frequency phenom- 
ena. KHz QPOs tend to have higher frequencies at higher countrates, but the correlation 
is not straightforward. On small timescales, in fact, the variation of the frequency with 
the countrate is usually very rapid, with the frequency going up of hundreds of Hz with 
slight increases of the countrate. If one plots the frequency of the oscillations versus the 
countrate, it is evident that the variation happens on different tracks, almost parallel to 
each other but distinct. Instead, there is a good correlation between the frequency and 
the spectral state [5], and at the lower and upper limits of their range the frequency dif- 
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FIGURE 1. In this figure it is possible to observe (left) the matter streams falling on the polar region 
from the funnel flow and the ones from instabilities at the inner radius and (right) the hotspots left on the 
surface by these two different accretion mechanisms [10] 



ference between the peaks decreases [6, e.g.]. Possible models for kHz QPOs include a 
number of different phenomena, including but not limited to relativistic resonances [e.g. 
7], Lense-Thirring relativistic precession of matter in preferred orbits [e.g. 8], magnetic 
field-driven oscillations [e.g. 9]. 

MOVING HOT SPOTS AND QPOS 

Global 3D magnetohydrodynamical simulations of accreting, weakly-magnetized neu- 
tron stars show (e.g. [1 1], [12], [13]) that most accretion configurations produce hotspots 
on the surface of the star. These spots are confined on a small region of the star (hence, 
by all means, "fixed") only in one case: when accretion is on the polar region through 
funnel flow and the misalignment angle between the magnetic field and the rotation axis 
is sufficiently large. In this case, in fact, matter finds a "preferred path" during its fall, 
that takes near the magnetic pole, fixed in the star's reference frame. This hotspot is thus 
tied to a region of the star and rotates with it, and an observer at infinity will see the 
light coming from the hotspot modulated at the rotational frequency of the star, simi- 
larly to what is observed in x-ray pulsars. In all other cases, namely funnel flow with 
small misalignment angle and unstable accretion in the equatorial region, the hotspots 
tend to move (Fig. 1): in fact, in these cases the magnetic field does not change con- 
siderably along and thus matter is able to fall in different places of the surface, or in 
different positions around the magnetic pole. If a stream of matter is formed that falls on 
the surface, the corresponding hotspot will move with an angular velocity similar to the 
velocity at the starting point of the stream in the disk. 

If we model with Monte Carlo techniques these moving hotspots, we extract a 
lightcurve and we take a Fourier Transform of it, we immediately find some interest- 
ing things. First of all, the shape of the peaks at the rotational frequency of the hotspot 
is a "bump", a QPO. This is because, even if the rotational frequency of the hotspots 
is fixed, their duration is limited and thus Fourier theory shows that the Fourier trans- 
form is not a single-frequency peak but a bell-shaped function (in the ideal case, a sine 
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FIGURE 2. Frequency variation of the hot spot movements with the mass accretion rate [10]. Results 
are scaled for the two cases B — 10 8 G and B = 5 • 10 8 G. The first case (right) shows only stable accretion 
from the funnel flow. The hot spot moves around the magnetic pole with a measurable angular velocity, 
which we plot against the mass accretion rate M. The second case (left) shows instabilities at the inner 
radius. In this case we have both accretion onto the magnetic pole with the usual funnel flow, like in the 
first case, and accretion from instabilities, moving approximately at the frequency of the inner disk. Hence 
the "double track" in the plot. 



function). Moreover, the rotational frequency of the hotspots is not perfectly stable, and 
thus the broadness of the peak increases. Second, the possible frequencies to observe at 
a time are two: one given by the hotspots moving around the magnetic pole, one given 
by the hotspots produced through instabilities. The frequency of the instability hotspots 
is that of the inner edge of the disk, not very stable, and the instabilities have a relatively 
short duration. The polar region hotspots, instead, are more stable and last longer. There- 
fore, the features they produce in power spectra are generally more coherent than those 
from instabilities. The frequency is lower because the accretion on the polar region does 
not come from the very edge of the disk, but from a larger region behind it. Finally, the 
moving hotspots produce red noise, whose shape and amplitude depend on the statistical 
distribution of amplitude and intensity of the hotspots. The striking similarities between 
the behavior of these two different ways to produce hotspots and the two kHz QPOs lead 
us to propose this mechanism as a possible explanation for the kHz QPO phenomenon. 



TESTING THE RELIABILITY OF THE MODEL 

To propose this mechanism as a viable way to explain QPO, we needed to be sure that 
the frequency of rotation of the hotspots followed a consistent behavior. First of all, the 
hotspots had to show up in different simulations with similar parameters in a similar 
way. Then, their frequencies had to change in a predictable way. 

For this reason, we performed a number of simulations differing only by the mass 
accretion rate. All the systems reproduced in the simulations have M = 1.4M , R = 
10km and misaligment angle G = 2°. One system has a rotational period T = 4.1ms, 
the other % — 5.4ms. In Fig. 2 we show the frequencies of all hotspots appearing on 
the surface, plotted against the mass accretion rate M. In one case the accretion is only 
through the funnel flow, in the other both through the funnel flow and the instabilities. 



As we can see from the figure, the frequencies of the hotspots in the funnel flow and 
from the instabilities increase with frequency, forming two distinct tracks in the v — M 
plane. Their frequency difference is around 200 Hz. Their consistent behavior tells us 
that these phenomena are something intrinsic to the accretion process as modeled in 
MHD simulations, and most likely to the accretion in real systems. 

CONCLUSIONS 

The results here reported represent a viable explanation for the kHz QPO phenomenon. 
To summarize the main points: 1) MHD simulations show that multiple, moving hotspots 
are formed on the surface of neutron stars, both in the funnel flow (if the misalignment 
angle is small enough) and from instabilities; 2) the movement of these hotspots pro- 
duces a modulation of the lightcurve, whose appearance on a power spectrum is a QPO 
plus red noise components; 3) the two different accretion mechanisms produce qualita- 
tively different QPOs, the ones from the funnel flow being usually more coherent and at 
a lower frequency; 4) their frequency separation is consistent with the frequency sepa- 
ration of real-life kHz QPOs. For more information, see [10]. 

ACKNOWLEDGMENTS 

Research supported by INFN (MB) and the contract ASI-INAF 1/088/06/0 for High En- 
ergy Astrophysics (MB). MMR was supported in part by NASA grant NNX08AH25G 
and by NSF grant AST-0807129. Simulations performed on NASA (Pleiades, 
Columbia), and CINECA (BCX, SP6, under the 2008-2010 INAF/CINECA agree- 
ment) supercomputers. 

REFERENCES 

1. T. Belloni, D. Psaltis, and M. van der Klis, ApJ 572, 392 (2002). 

2. M. van der Klis, J. H. Swank, W. Zhang, K. Jahoda, E. Morgan, W. H. G. Lewin, B. Vaughan, and 
J. van Paradijis, ApJ 469, LI (1996). 

3. T. di Salvo, M. Mendez, and M. van der Klis, A&A 406, 177 (2003). 

4. D. Barret, J.-F. Olive, and M. C. Miller, MNRAS 370, 1 140 (2006). 

5. M. Mendez, M. van der Klis, E. C. Ford, R. A. D. Wrjnands, and J. van Paradijis, ApJ 511, L49 
(1999). 

6. S. Boutloukos, M. van der Klis, D. Altamirano, M. Klein-Wolt, R. A. D. Wrjnands, P. G. Jonker, and 
R. P. Fender, ApJ 653, 1435 (2006). 

7. W. Kluzniak, and M. A. Abramowicz, eprint arXiv p. 5057 (2001). 

8. L. Stella, and M. Vietri, ApJL 492, L59 (1998). 

9. L.-X. Li, and R. Narayan, ApJ 601, 414 (2004). 

10. M. Bachetti, M. M. Romanova, A. Kulkarni, L. Burderi, and T. di Salvo, MNRAS 403, 1 193 (2010). 

11. M. M. Romanova, G. V. Ustyugova, A. V. Koldoba, and R. V. E. Lovelace, ApJ 610, 920 (2004). 

12. A. K. Kulkarni, and M. M. Romanova, MNRAS 398, 701 (2009). 

13. M. M. Romanova, and A. K. Kulkarni, MNRAS 398, 1 105 (2009). 

14. M. T. Menna, L. Burderi, L. Stella, N. Robba, and M. van der Klis, arXiv astro-ph (2003), 

astro-ph/0303606vl. 



